NASA-CR-204621 




/ S 




Effect of Correlated Precision Errors 
on Uncertainty of a Subsonic Venturi Calibration 

S. T. Hudson, W. J. Bordelon Jr., H. W. Coleman 


Reprinted from 

DIMM 

Volume 34, Number 9, Pages 1862-1867 


&AIAA 

A publication of the 

American Institute of Aeronautics and Astronautics, Inc. 
1801 Alexander Bell Drive, Suite 500 
Reston.VA 22091 



AIAA Journal 

Vol. 34, No. 9, September 1996 


Effect of Correlated Precision Errors on Uncertainty 
of a Subsonic Venturi Calibration 

S. T. Hudson* and W. J. Bordelon Jr. f 
NASA Marshall Space Flight Center, Huntsville, Alabama 35812 

and 

H. W. Coleman* 

University of Alabama in Huntsville , Huntsville, Alabama 35899 


An uncertainty analysis performed in conjunction with the calibration of a subsonic venturi for use in a turbine 
test facility produced some unanticipated results that may have a significant impact in a variety of test situations. 
Precision uncertainty estimates using the preferred propagation techniques in the applicable American National 
Standards Institute/American Society of Mechanical Engineers standards were an order of magnitude larger than 
precision uncertainty estimates calculated directly from a sample of results (discharge coefficient) obtained at the 
same experimental set point The differences were attributable to the effect of correlated precision errors, which 
previously have been considered negligible. An analysis explaining this phenomenon is presented. The article is not 
meant to document the venturi calibration, but rather to give a real example of results where correlated precision 
terms are important The significance of the correlated precision terms could apply to many test situations. 


Nomenclature 

B = bias limit estimate 

C d - discharge coefficient 

d ~ diameter 

J = number of measurement variables 
M - number of tests (measurement sets) 
M n = Mach number 

N = number of readings 

P = precision limit estimate 

p = pressure 

Re = Reynolds number 

r = result 

S = sample standard deviation 

T = temperature 

U = uncertainty estimate 

tV = mass flow rate 

X = measurement reading 

p = bias component of error 

8 = total error 

e = precision component of error 

Q = sensitivity coefficient 

/i = population mean 

p - correlation coefficient 

a = population standard deviation 

Subscripts 

0 = total 

1 = venturi inlet 

2 = venturi throat 
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Introduction 

P RECISION uncertainty estimates using the preferred propaga- 
tion techniques in the applicable American National Standards 
Institute/ American Society of Mechanical Engineers (ASME) stan- 
dards have traditionally assumed that the correlated precision terms 
are negligible. This assumption has been made because precision 
errors are considered to be random. However, the randomness as- 
sumption means that there is no correlation among sequential mea- 
surements of a single variable. But the real question concerns the 
correlation of precision errors among simultaneous measurements 
of different variables. This article presents results from a subsonic 
venturi calibration as a real example of a case where the correlated 
precision terms are significant in the uncertainty analysis. 

Mass flow venturis are used in a variety of test situations to deter- 
mine the mass flow rate of gases or liquids. To accurately determine 
the mass flow rate using a venturi, the discharge coefficient must be 
known. The discharge coefficient is defined as the ratio of the actual 
mass flow passing through the venturi to the ideal one-dimensional 
inviscid mass flow rate. Critical flow venturis are often used because 
Cd is well known for standard designs such as the Smith-Matz 1 and 
the ASME ! 2 nozzles. However, a critical flow venturi is not al- 
ways the most desirable option. 2 When subsonic flowmeters must 
be used, determination of the discharge coefficient is crucial to meter 
accuracy. 

A subsonic Herschel-type venturi is used in the Turbine Test 
Equipment (TTE) at the Marshall Space Flight Center (MSFC). 3 
The subsonic venturi was chosen over a critical flow venturi to 
minimize flow noise. The TTE is used to test liquid rocket engine 
turbines in air at scaled conditions. A technology demonstrator gas 
generator oxidizer turbine currently is being tested in the TTE. 4 
The test data from this technology turbine are to be used for code 
validation as well as performance evaluation; therefore, absolute ac- 
curacy is very important. The mass-flow-rate measurement is used 
to calculate the turbine efficiency and the turbine flow parameter. 
The goal of the technology turbine testing is to determine these two 
performance parameters within ±1%. This means that an accuracy 
of ±0.5% is needed for the mass flow measurement. The accuracy of 
the mass flow measurement in the TTE depends on the accuracy of 
the venturi inlet pressure, throat pressure, and inlet temperature mea- 
sured during the test; the venturi inlet and throat diameter measured 
during inspection; and the venturi discharge coefficient determined 
during calibration. To meet the ±0.5% goal for the TTE mass flow 
measurement, the discharge coefficient from calibration needs to be 
determined within approximately ±0.25%. 
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The TTE venturi had been calibrated, but the stricter requirements 
for the technology turbine testing prompted a review of the original 
calibration. This review revealed several concerns. First, the tech- 
nology turbine would be operating outside the calibration range for 
the majority of the test matrix. Also, the uncertainty of the standard 
nozzle used to determine the calibration mass flow rate was ±0.5%. 
The dischaige coefficient uncertainty would thus have been greater 
than the required ±0.25%. 

On the basis of the above findings, it was decided that a new cali- 
bration of the TTE venturi was needed. The question was whether or 
not the venturi could be calibrated to the required ±0.25%. A nine- 
point test matrix consisting of three venturi throat Mach numbers 
(0.2, 0.5, and 0.7) and three venturi throat Reynolds numbers (10 6 , 
3 x 10 6 , and 6 x 10 6 ) was decided upon. FluiDyne Aerotest Group 
of Aero Systems Engineering, Inc. (ASE), was asked to evaluate 
the uncertainty in C d that they could obtain in their test facility for 
this test matrix. ASE estimated a bias error and a precision error 
for each measurement at each test point. These estimates were then 
propagated through the data reduction equation for C d to estimate 
an uncertainty in the discharge coefficient. On the basis of these 
estimates, it appeared that the calibration goal of ±0.25% could be 
achieved at Mach numbers of 0.5 and 0.7 with very careful calibra- 
tion. However, the low Mach number (0.2) points seemed to be a 
problem. Large precision errors were predicted at these points be- 
cause of the low pressure differential between the venturi inlet and 
throat at these conditions. On the basis of past experience, however, 
ASE felt that they could achieve better precision with their test sys- 
tem than that predicted. It was decided to proceed with the calibra- 
tion and to determine the precision limit in the discharge coefficient 
directly from the test data. Precision limits determined from the test 
data then would be combined with the bias estimates to determine an 
overall uncertainty for each test point. These uncertainties would be 
compared with the pretest predictions and calibration requirement. 

When the test data were obtained, two methods of estimating the 
precision limits in the calibration data were used: 1) calculating the 
standard deviations for the pressure and temperature measurements 
and propagating these through the data reduction equation while 
assuming uncorrelated precision errors and 2) calculating the stan- 
dard deviations directly from multiple C d results. The actual test 
data were used in both cases so that a direct comparison could be 
made. Method 2 gave precision limits that were an order of mag- 
nitude less than those calculated using method 1 . The calibration 


approach, calibration results, and an uncertainty analysis explain- 
ing the differences between the two methods of determining the 
precision limits are discussed. 

Calibration Approach 

The venturi calibration runs were conducted in air in the Chan- 
nel 7 test facility at ASE Medicine Lake Aerodynamic Laboratory. 
The calibration test runs were conducted between Jan. 2 1 and Feb. 1 , 
1994. This section describes the experimental setup, measurements, 
and test procedure as well as details of the test venturi itself. More 
detailed information is presented in the calibration report. 5 

The venturi used in the MSFC TTE is a Herschel type as defined 
by Bean. 6 The nominal inlet diameter is 7.980 in., and the nominal 
throat diameter is 3.396 in. This venturi was sized to operate at 
subsonic conditions to minimize flow noise or pressure fluctuations 
caused by a normal shock. The venturi is made of stainless steel 
and has two instrumentation planes. The venturi inlet plane contains 
two gas temperature probes and four static pressure taps. The venturi 
throat contains four static pressure measurements. The overall length 
of the venturi is 59 in. 

The facility used for the venturi calibration was Channel 7 at ASE 
modified as shown in Fig. 1. The facility consists of a 500-psi air 
storage source, heater, inlet control valve, ASME metering nozzle, 
settling chamber, test venturi, and backpressure control valve. The 
air supply, heater, control valve, metering nozzle, and settling cham- 
ber are part of the existing facility setup. An adapter was fabricated 
to attach the venturi to the settling chamber discharge. Downstream 
of the venturi, a 16-6-in. reducer was made from standard pipe re- 
ducers to install the 6-in. backpressure control valve. The flow was 
exhausted to atmospheric conditions downstream of the backpres- 
sure control valve. Two ASME long-radius metering nozzles were 
used as the calibration standards to cover the range of calibration 
conditions. The ASME nozzle nominal throat diameters were 1 .47 
and 1.75 in. Both nozzles were operated at choked conditions for 
all test runs. 

Pressure and gas temperature measurements were recorded for 
both the metering nozzle (calibration standard) and the test venturi. 
The inlet total pressure, two inlet gas temperatures, and a static pres- 
sure just downstream of the throat were measured on the metering 
nozzle. Measurements taken on the test venturi included four inlet 
static pressures, four throat static pressures, and two inlet gas tem- 
peratures. Atmospheric pressure also was measured for each run. 



Fig. 1 Venturi calibration schematic (not to scale). 
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The pressure measurements were made with a Pressure Systems, 
Inc., pressure scanning system. The metering nozzle gas tempera- 
tures were measured with shielded type J thermocouples, and the 
venturi gas temperatures were measured with type K thermocou- 
ples. The voltage outputs from the thermocouples were measured 
using an electronic digital data acquisition system. 

The venturi was calibrated over the range of throat Mach num- 
bers and throat Reynolds numbers by setting the metering nozzle 
upstream pressure (mass flow for choke condition) with the upstream 
control valve and the test venturi backpressure with the backpressure 
control valve. With the metering nozzle choked, the test mass flow 
and venturi Reynolds number were fixed for a constant metering 
nozzle inlet total pressure and temperature. The venturi backpres- 
sure then was adjusted to set the venturi Mach number. All test 
runs were conducted with the air at room temperature. Once the test 
conditions were established, the pressure and temperature measure- 
ments were scanned 10 times with about 6 s between each scan. 


95% confidence estimate of P is interpreted to mean that the ±P 
interval about a single reading of X ) should cover the (biased) parent 
population mean /i 95 times out of 100. 

In nearly all experiments, the measured values of different vari- 
ables are combined using a data reduction equation (DRE) to form 
some desired result. A general representation of a DRE is 

r = r(X {9 X 2 ,...,Xj) (1) 

where r is the experimental result determined from J measured 
variables X, . Each of the measured variables contains bias errors and 
precision errors. These errors in the measured values then propagate 
through the data reduction equation, thereby generating the bias and 
precision errors in the experimental result r. 

If the large sample assumption is made, 8 then the 95% confidence 
expression for U r becomes 

Uf = B; + p} (2) 


Calibration Results 

The venturi was calibrated over a nine-point test matrix. The 
nominal throat Mach numbers for this matrix were 0.2, 0.5, and 
0.7. The nominal throat Reynolds numbers for this matrix were 
10 6 , 3 x 10 6 , and 6 x 10 6 . Ten data scans were taken at each set 
point. The venturi pressures and temperatures were recorded and 
used to calculate Mn, Re, and C d for each scan. The equations for 
these calculations are closed form and are given in the Appendix. 6 7 
The 10 scans then were averaged to obtain mean values of Mn , Re, 
and C rf . Three runs were made (10 scans each, except run 808, which 
had 9 scans) for the worst-case data point {Mn — 0.2, Re = 10 6 ). 
This was the point where the uncertainty was expected to be the 
highest. These repeat runs were made to provide more data for 
understanding the precision behavior associated with this point. The 
average values of Mn, Re, and C d for each run are given in Table 1 . 
(Note: At Mn = 0.7, Re had to be increased to 1 .5 x 1 0 6 to achieve 
the test point.) 


Uncertainty Analysis Overview 

The word accuracy is generally used to indicate the relative close- 
ness of agreement between an experimentally determined value of 
a quantity and its true value. Error (<5) is the difference between 
the experimentally determined value and the truth; thus as error de- 
creases, accuracy is said to increase. Only in rare instances is the 
true value of a quantity known. Thus, it is necessary to make an 
estimate of an error, and that estimate is called an uncertainty, U . 
Uncertainty estimates are made at some confidence level. A 95% 
confidence estimate, for example, means that the true value of the 
quantity is expected to be within the ±U interval about the experi- 
mentally determined value 95 times out of 100. 

Total error 8 can be considered to be composed of two compo- 
nents: a precision (random) component e, and a bias (systematic) 
component p. An error is classified as precision if it contributes 
to the scatter of the data; otherwise, it is a bias error. As an esti- 
mator of p, a systematic uncertainty or bias limit B is defined. A 
95% confidence estimate is interpreted as the experimenter being 
95% confident that the true value of the bias error, if known, would 
fall within ±B. As an estimator of the magnitude of the precision 
errors, a precision uncertainty or precision limit P is defined. A 


where we define the bias limit (systematic uncertainty) of the result 
as 


j j- i j 

B r = Y 9 ~ B ‘ +2 Y Y BAp^B.B, (3) 

1-1 »' = ] k =1 + 1 

and the precision limit (precision uncertainty) of the result as 

j J~\ j 

P r = Y 0 ? P ? + 2 Y Y W 

; = 1 | = 1 k « I + 1 


where pbik is the correlation coefficient appropriate for the bias errors 
in X t and X k , p sUc is the correlation coefficient appropriate for the 
precision errors in X { and X k , and 


e { = 


dr 

dXj 


(5) 


The 95% confidence precision limit for a variable X, can be esti- 
mated as 


Pi = 2 Si 

where the sample standard deviation for X { is 

i , . i 

and the mean value for X, is defined as 


( 6 ) 


(7) 


X 


1 

77 


N 




( 8 ) 


Traditionally, the correlation terms in Eq. (4) have been assumed 
to be negligible. This assumption has been made because precision 
errors are considered random. However, the randomness assump- 
tion means that there is no correlation in sequential measurements 
of a single variable. But, the real question concerns the correla- 
tion of precision errors in simultaneous measurements of different 
variables. These correlation terms may be significant after all. 


Table 1 Calibration results 


Run no. 

Mn 

Re x 1CT 6 

Q 

1.01-1.10 

0.197 

0.975 

0.9899 

1.11-1.20 

0.192 

0.966 

0.9917 

808.01-808 10 

0.201 

1.084 

0.9866 

4.1 1-4.20 

0.199 

2.932 

0.9894 

7.01-7.10 

0.199 

5.953 

0.9933 

2.06-2.15 

0.502 

0.977 

0.9885 

5.01-5.10 

0.495 

2.971 

0.9914 

8.01-8.10 

0.493 

5.834 

0.9927 

3.11-3.20 

0.698 

1.450 

0.9907 

6.1 1-6.20 

6.695 

2.951 

0.9914 

9.01-9.10 

0.683 

5.894 

0.9938 


Single Test: Single Readings 

The precision limits defined in Eq. (6) and used in Eq. (4) are 
applicable to a single test with single readings — that is, at a given 
test condition, the result is determined once using Eq. (1), and the 
X i are considered single measurements. This is the situation often 
encountered in large-scale engineering tests in which measurements 
of the variables are made at a given set point over a period that is 
small compared to the periods of the factors causing variability in 
the experiment. A proper precision limit (one indicative of the dis- 
persion of the variable over several cycles of the factors causing 
its variation) cannot be calculated from readings taken over such 
a small time interval. For such data, the measurement(s) of a vari- 
able Xi should be considered a single reading — whether the value 
of X , is the average of 10, 10 3 , or 10 6 readings taken during the 
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short measurement time. In such a test, the value for the precision 
limit to be associated with a single reading would have to be based 
on previous information about that measurement obtained over the 
appropriate time interval. 9 If previous readings of a variable over an 
appropriate interval are not available, then the experimenter must 
estimate a value for P, using the best information available at that 
time. 8 * 10 

Multiple Tests: Single Readings 

If a test is performed so that M multiple sets of measurements 
(X| , X 2 , . . . , Xj)k at the same test condition are obtained, then M 
results can be determined using Eq. (1), and the best estimate of the 
result r would be r , where 


r = 



(9) 


If the M (M > 10) sets of measurements were obtained over an 
appropriate time period, the precision limit that should be associated 
with this averaged result would be 

P- r = 2S r /-/M (10) 


and for a single result would be 

P r — 2S r (11) 


where S r is the standard deviation of the sample of M results 


S r = 


— E< r ‘ - ?)2 


M - 1 


( 12 ) 


Single Test: Averaged Readings 

Some previous approaches 1 1-13 considered the general case to 
be that in which the best estimate of the result is determined using 
averages of the X, so that 

r = r(X, t X 2 Xj) (13) 

and specified that (for large samples) the precision limits be taken 
as P,/(Ni ) l/2 in Eq. (4) (with the correlation terms neglected). The 
(p sik P, P k ) terms in Eq. (4) take into account the possibility of pre- 
cision errors in different variables being correlated, and these terms 
have traditionally been neglected, as mentioned previously. How- 
ever, precision errors in different variables caused by the same un- 
controlled factor(s) are certainly possible, as shown in this article. 
In such cases, one would need to acquire sufficient data to allow a 
valid statistical estimate of the correlation coefficients to be made 
in Eq. (4). Note, however, that the approach using Eqs. (10-12) im- 
plicitly includes the correlated error effect — a substantial advantage 
over the approach of Eq. (1 3). 

Uncertainty Analysis Results 

The uncertainty analysis presented here focuses on precision un- 
certainties. Two methods were used to evaluate the precision limit 
associated with the experimental result, C d . First, the precision limit 
of the result was calculated by propagating the precision of the av- 
eraged measured variables through the data reduction equation, as 
discussed in the paragraph above, while neglecting correlated preci- 
sion errors (method 1). Equations (4-8) and (13) were used for this 
approach. Second, the precision limit of the result was calculated 
directly from the sample standard deviation of C d for each test point 
(method 2). Equations (9-12) and (1) were used for this approach. 
Again, actual test data were used in both cases so that a direct com- 
parison between the two methods could be made. The two methods 
are discussed, and the results are compared and explained in the 
following paragraphs. 

A Fortran program was written to calculate the precision limit of 
C d by the propagation method (method 1 ). It was assumed that only 
/?, , p 2j and T\ contributed to the precision uncertainty: 

Pc„ = Pc d (Pu Pi,T\) 0 4 ) 


Table 2 C d precision results 


Mn 

Re x 10“ 6 

Cd 

P C,.% 

Method 1 Method 2 

0.197 

0.975 

0.9899 

2.00 

0.12 

0.192 

0.966 

0.9917 

4.37 

0.1 1 

0.201 

1.084 

0.9866 

2.15 

0.043 

0.199 

2.932 

0.9894 

2.46 

0.082 

0.199 

5.953 

0.9933 

1.17 

0062 

0.502 

0.977 

0.9885 

0.37 

0.050 

0.495 

2.971 

0.9914 

0.75 

0.053 

0.493 

5.834 

0.9927 

0.14 

0 061 

0.698 

1.450 

0.9907 

0.09 

0.047 

0.695 

2.951 

0.9914 

0.23 

0.028 

0.683 

5.894 

0.9938 

0.07 

0.024 


The precision uncertainties associated with the mass flow rate deter- 
mined by the standard and the test venturi inlet and throat diameters 
were assumed to be negligible relative to the precision uncertainties 
of the test venturi measurements of pressure and temperature. The 
standard deviations 5 of the three measurements ( p\ , P 2 , and TO 
were calculated from the 10 data scans for each test run (N = 10). 
The precision limit for each of these mean quantities was then cal- 
culated using P = 25/(10) 1/2 . The values of p u Pi. and T\ were 
set to the average of the 10 scans for each point. The code used these 
average values and the precision estimates as inputs for the variables 
and propagated these values through the DRE to obtain an estimate 
of the precision limit for C d . A central-difference numerical scheme 
was used to approximate the partial derivatives. No correlation of 
the precision errors in the measured variables was considered. 

A summary of the results using the propagation method 
(method 1) is given in Table 2. These results indicate that the goal 
of achieving an uncertainty in C d of ±0.25% cannot be attained. 
Notice particularly the extremely high precision uncertainties cal- 
culated for the low Mach number points. 

The precision limit of the discharge coefficient was then deter- 
mined directly from the C d results (method 2). C d was calculated for 
each data scan. A sample standard deviation 5 was then computed 
for each test point using the 10 data scans (M = 10). The precision 
was again calculated using P - 25/(10) 1/2 . These results also are 
given in Table 2. The results indicate that the goal of an uncertainty 
in C d of ±0.25% can be attained. 

The precision limits calculated by method 1 were much higher 
than those calculated by method 2. The differences were especially 
large at the low Mach number points. The large uncertainties at 
the low Mach numbers were attributable to the small difference be- 
tween the venturi inlet and throat pressures at these conditions. An 
examination of the influence coefficients obtained from the propa- 
gation method revealed that the C d precision limit was extremely 
sensitive to the venturi inlet and throat pressure measurements. This 
sensitivity to the static pressure is discussed elsewhere. 2 The venturi 
temperature measurement had very little effect on the precision of 
Crf. The venturi mass flow and, therefore, C d are functions of the 
ratio of the inlet to the throat pressure. This means that variations in 
the inlet pressure and the throat pressure will not affect the precision 
of W or C d if the ratio remains constant. 

The propagation method treated the precision errors in the two 
venturi pressures as independent. A plot of these two pressures nor- 
malized to the critical flowmeter inlet total pressure for a particular 
run (Fig. 2) shows that the variations of the two pressures are not 
independent, however. This same trend was seen for all test condi- 
tions. The fact that the throat pressure varied with the inlet pressure 
was a function of the test facility control. The distance between 
the critical flowmeter and the venturi was small; therefore, pres- 
sure variations in the critical flowmeter and in the venturi were in 
phase. The variations in the pressure measurements were not truly 
random; they were correlated. This correlation was not accounted 
for in the propagation analysis. Method 2 automatically accounted 
for the correlation effect because C d was calculated for each scan, 
and these data were used to calculate the precision limit. To ac- 
count for the correlation effect in the propagation analysis, a corre- 
lation coefficient must be estimated. To use this method for pretest 
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Fig. 2 Normalized venturi inlet and throat pressures. 


Table 3 Multiple test results 


Run 

Mn 

Re x l<r 6 

Q 

Pc d .% 

1.01-1.10 

0.197 

0.975 

0.9899 

0.12 

1.11-1.20 

0.192 

0.966 

0.9917 

0.11 

808.01-808.10 

0.201 

1.084 

09866 

0.043 

Overall 

0.197 

1.008 

0.9895 

0.10 


predictions, one must recognize that the precision errors of the two 
pressures will be correlated and include the correlation terms in the 
analysis. 

Note that the precision limit should be calculated using both meth- 
ods. Often, the precision limit using method 1 is computed automat- 
ically in the data reduction program. When this is done, it is a fairly 
simple task to add the method 2 calculation of the precision limit to 
the program. Having the precision limits calculated by both methods 
would allow one to make a comparison between the two methods. If 
the precision limits from the two methods are significantly different, 
then the possibility of having correlated precision effects should be 
evaluated. 

To gain more insight into the precision of the worst-case point of 
Mn = 0.2 and Re = 10 6 , three runs were made at this test condition. 
Two of the runs (10 scans each) were made at different times on 
the same day, and the third run (9 scans) was made almost two 
weeks later. The data for each run were calculated by method 2, as 
discussed previously. Combining the three runs gave M = 29 for 
the overall data set. These data are summarized in Table 3. The data 
show that the run-to-run repeatability is good. Some variation in C d 
is expected because of set-point variations. The precision limit of 
the overall average discharge coefficient is 0.10% when all three 
runs are considered [P = 2S/(29) 1/2 ]. This data set allows one to 
be more sure of the precision of the worst-case point. 

Conclusions 

Correlated precision terms can have a significant effect on the pre- 
cision limit of a result. Precision uncertainty estimates using prop- 
agation techniques without considering correlated precision terms 
were an order of magnitude larger than precision uncertainty es- 
timates calculated directly from a sample of results obtained at 


the same experimental set point. The potential effect of correlated 
precision terms makes it better to use multiple results rather than the 
propagation method to determine a precision limit for a result. The 
propagation method is useful in the planning phase, but one must 
recognize the effect of the correlated precision terms and include 
this in the analysis if it is applicable in a particular test situation. 
When the precision limit is automatically computed in the data re- 
duction program, it should be computed by both methods so that 
a comparison can be made. If the precision limits are significantly 
different, then the possibility of having correlated precision effects 
should be considered. The correlated precision effect is probably 
much more prevalent than considered in the past, because it occurs 
in any system in which multiple variables drift in unison in response 
to unsteadiness during the testing period. 
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Appendix: Equations 

The equations used to calculate the venturi Mach number, 
Reynolds number, mass flow rate, and discharge coefficient for this 
application are included here. The basic equations came from Bean, 6 
and they are applied to this specific case elsewhere. 7 (Note that the 
equations for Mach number and Reynolds number are the same as 
those used by ASE, 5 but the equation for mass flow is different.) The 
equations are written for English units; therefore, pressure must be 
in psia, area in square inches, and temperature in °R. Constants used 
include y = 1 .4 and R = 53.35 ft-lb//lb m °R. The symbol p is now 
used for density ( p was the correlation coefficient in the text of the 
article). Also, the symbol M (not Mn) \% now used for Mach num- 
ber. The constant 9.9 x 1 0 6 in Eq. (A7) is the thermal expansion 
coefficient. 
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Venturi diameter ratio and pressure ratio: 

DR = d 2 /di ; PR^Pi/Pi 


Throat Mach number: 


Af 2 = 


PICj_ - 1 
1 - ( DR*)(PR “) 


Inlet Mach number: 

Mi = DR 2 PR* Mi 


(Al) 


(A2) 


(A3) 


Inlet static temperature: 

T\ = 

Inlet density: 


1 4- 0.2M? 


Pi = 


144 pi 

^rT 


(A4) 


(A5) 


Throat adiabatic wall temperature: 


T _ T / 1 + 01784 Mj \ 

T2aw ~ To ' V 1 + 0 . 2 M| J 


Throat diameter corrected because of difference between operat- 
ing temperature and room temperature: 

d 2r = d 2 [ 1 + 9.9 x 10- 6 (r ju , - 529.67)] (A7) 

Throat area based on throat diameter in Eq. (A7): 

A 2t =7t(d 2T /2) 2 (A8) 

Throat Reynolds number per foot per psia: 

Re _ (181,176,192)M 2 [r o , + 198.6(1 + 0.2M?)] 
dpoi~ (r„,) 2 (i + 0 . 2 M 2 2 ) 25 


Throat Reynolds number: 

Re = Re 2 = ( Re/dp fn )p 2 (\ + Q.2Mlf' 5 (d 2T /\2) (A10) 


Ideal mass flow rate (lb m /s): 

= 1 250605(A 2r ) ^ t _ ^ 


(All) 


Discharge coefficient: 

= Wsundjri = Jsi ( A12) 

w idcal 
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